Pairing and shell evolution in neutron rich nuclei 
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From the experimental data on odd-even staggering of masses, we have shown that variation of 
pairing as a function of neutron number plays an important role in many distinctive features like 
occurrence of new shell closures, shell erosion, anomalous reduction of the energy of the first 2 + 
state and slower increase in the B(E2, 2f — >• Of) in the neutron-rich even-even nuclei in different 
mass regions. New predictions have been made in a model independent way. Results of theoretical 
calculations support the phenomenological findings. 

PACS numbers: 21.60.Cs,21.30.Fe,23.20.Lv,27.60.+j 



The energy of the first 2 + state (E(2f)) of even-even 
nucleus is a sensitive probe to study the evolution of shell 
structure. For spherical nuclei, this state is formed by 
breaking a pair spending the pair binding energy, show- 
ing pronounced maxima at the shell closures. Nuclei 
away from the shell closures are gradually driven away 
from spherical symmetry, where "anomalously" low first 
excited collective 2 + states are observed. 

Usually, a few distinctive features have been identified 
as the signatures of changing shell structure away from 
stability Three new doubly magic Oxygen isotopes 
have been observed. For neutron rich nuclei, semi-magic 
ones like 32 ,A4g20, io^e2o have shown erosion of N=20 
shell closure with sudden decrease in their E(2f) values 
and an increase in their corresponding B(£72, 2^ — >• Of ). 
However some recent measurements of reduced (E(2f)) 
and B(E2) values in highly neutron nuclei having Z=4-18 
have raised a serious discussion on hindered E2 strength 
unexpected for these nuclei phenomenologically as well as 
theoretically . In heavier nuclei, after observation Q 
of reduction of both E(2+) and B(E2) in 136 Te, Terasaki 
et al. [7], from QRPA calculations, traced the origin of 
this anomalous behaviour in 136 Te isotope to a reduced 
neutron pairing above the N=82 magic gap. Around this 
period, new empirical interactions [8j proposed for neu- 
tron rich isotopes above the 132 Sn core and later those 
for neutron rich nuclei in the sd — fp shell [t| [l(| also 
included reduction of pairing matrix elements for better 
reproduction of data. So these studies indicate an im- 
portant role of pairing in the evolution of the structure 
of exotic nuclei. 

In this letter, we show from the neutron and proton 
pairing energies estimated from the differences in exper- 
imental binding energies in different mass regions, that 
pairing plays a crucial role in evolution of shell structure 
in atomic nuclei. For observed shell erosions at N=8 in 
Be isotopes and at N=20 in Ne, Mg isotopes, relative 
enhancement of proton pairing along with weakening of 
neutron pairing have been found to be important. We 
have found a clear correlation among the occurrence of 
new shell closures, the reduction of £7(2*) and hindered 



FIG. 1: Variation of (a) A n for odd isotopes of Be, A p for 
odd isotopes of Li and B, (b) E(2f) energies for even isotopes 
of Be with increasing neutron number (N). 



B(E2) in the neutron-rich even-even nuclei and weaken- 
ing of neutron pairing. This feature has been corrobo- 
rated by Shell Model (SM) and Cranked Hartree Fock 
Bogoliubov (CHFB) results. 

The empirical neutron and proton pairgaps are related 
to the odd-even staggering of binding energies [ll| . They 
are defined as 

A n (Z,N) = ^-[B(Z,N-l)-2*B(Z,N) + B(Z,N + l)} 
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[B(Z - 1, N) - 2 * B(Z, N) + B(Z + 1, N)] 



where B(Z,N) is the binding energy [12j of the nucleus 
with Z protons and N neutrons. The factor depending on 
the number parity 7Tjv (tz) is chosen so that the pairing 
gaps are all positive. To exclude the effect of variation 
of single particle energies in the staggering, the pairing 
gaps for odd N (Z)'s have been considered as the neutron 
(proton) pairing energies. The neutron (proton) pairing 
energy at an even N (Z) is estimated from the average 
of corresponding values for two neighbouring odd N (Z) 
nuclei. 

In Be (Z=4) isotopes, the N=8 shell gap quenching 
has been observed [J]. Comparison of Figs. [TJt and b 
reveals how 2 + states are built in Be isotopes as N in- 
creases. In 6 Be, this state is built by breaking a weak 
7rlp 3 /2 proton pair. In s,la Be, both protons and neutrons 
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FIG. 2: Variation of (a) A n for odd isotopes, (b) E(2f) en- 
ergies for even isotopes, of C and O with increasing neutron 
number (N). 



can contribute. The 2j state in 10 Be can be generated 
by breaking a neutron pair in and exciting one 

of them to nlpi/ 2 . As a result the 2+ energy increases 
(Fig. [Tp). With increasing neutron number, the neu- 
tron pairing decreases with increase in proton pairing. 
At around N=8, the proton pairing becomes stronger 
(Fig. QJi) by about 2 MeV. In 12 Be, instead of breaking 
a strong proton pair (A p ~ 3.1 MeV) in 7rlp 3 / 2 or excit- 
ing two neutrons across the N=8 shell closure to vld 5 / 2 , 
the most favourable option to generate a 2 + state, is to 
deform the system (nuclear Jahn Teller effect [l]|) and 
bring down the energy of the vld 5 / 2 orbit. The orbit 
below the gap, i.e the 1/2[101] Nilsson orbit originating 
primarily from vlpi/ 2 at low deformations, is insensitive 
to increase in deformation compared to the downslope of 
the 1/2 [220] orbit of vld b/2 [14]. This results in closing 
of the N=8 shell gap for 12 Be isotope. The single particle 
energy (spe) difference v(lpi/ 2 ~ 1^5/2) m Be, at N=8, 
estimated from difference of A„ of relevant even N and 
(N+l) [11] isotopes is w 1.8 MeV compared to » 3.9 MeV 
in corresponding C isotopes providing evidence of shell 
erosion. The weakening of the neutron pairing (A„ ~ 0.8 
MeV at N=10) further decreases the E(2f) energies of 
14 Be without substantial increase in deformation j4|. 

For C and O isotopes at N=8, the proton pairing is 
nearly equal to the neutron pairing. Fig. [2^, shows that 
the A„'s for carbon and oxygen isotopes decrease with 
N above N=ll. For the closure of the neutron 1^5/2 
and 2si/2 subshells at N=14 and 16, the semi-magic oxy- 
gen isotopes exhibit shell closures at 22 ' 24 0, respectively 
(Fig. [2p). The presence of two valence proton holes in- 
hibits the occurrence of shell closure for the carbon iso- 
topes. The spe difference v{\d^i 2 — 2s 1 / 2 ) and A„ at 
N=14, are ~ (2.7, 1.1) MeV and (1.3, 1.0) MeV for O and 
C isotopes, respectively. In 24 0, v{2s\/ 2 ~ 1^3/2) is — 3.9 
MeV compared to A„ ~ 0.6 MeV. However, the weaken- 
ing of pairing in carbon isotopes is manifested through 
decrease in its E (2^) values. 

Another prominent example of the erosion of magicity 
is observed at N = 20 for Z = 10, 12. In Fig. rja, vari- 
ation of A p - A„ (=A), the difference between proton 
and neutron pairing energies with increasing N has been 
shown. For Ne and Mg isotopes, across N=20, the proton 
pairing shows an increase with corresponding decrease in 
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FIG. 3: Variation of (a) A p - A n and (b) E(2f) energies 
for even isotopes of TVe, Mg, Si, S and Ar with increasing 
neutron number (N). 
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FIG. 4: Same as Fig. [2] for Sn and Te. "New" indicates the 
An calculated from most recent mass measurements Hal. 



neutron pairing (A positive). The situation is similar 
to Be isotopes at N=8. The Nilsson orbitals involved 
are upsloping 3/2 [202] from vld 3 / 2 and relatively faster 
downsloping 1/2 [330] from ^l/V/2- So onset of deforma- 
tion along with reduced neutron pairing depresses the 
E(2j) in N=20 isotopes of Ne and Mg without substan- 
tial increase in B(E2)s. The weakening of the neutron 
pairing continues for higher N with further decrease in 
the E(2f) energies of 32 iVe, 34 Mg. For Si, S and Ar nu- 
clei, below N=28, A is negative (except for Ar at N=16), 
indicating stronger neutron pairing. However, in Si nu- 
clei, A crosses the zero line and becomes positive from 
N~ 27, while in S and Ar crossing is at N=28. So Si shows 
an erosion of the N=28 shell closure manifested through 
depression in E(2^" ) energy (Fig.[3b) due to onset of de- 
formation where neutrons from upsloping 7/2[303] from 
vlf-j/2 can be excited to downsloping 1/2[321] of v2p 3 / 2 , 
with 5/2[312] from fl/7/2 being insensitive to increase in 
deformation 3]. In S and Ar isotopes, erosion is not dis- 
tinct, but beyond N=28, depression in E(2^~) with slower 
increase m B(E2) is observed @. 

The Sn region is even more interesting, showing how 
new features develop as one moves far away from sta- 
bilty for heavy nuclei. For neutron-rich isotopes of Sn 
and other heavier nuclei, neutron pairing shows reduc- 
tion. However, the proton pairing is generally smaller 
than or close to neutron pairing for these nuclei. The 
region above neutron-rich doubly closed 132 Sn has great 
similarities [15] with that above doubly closed 208 Pb on 
the stability. The even Sn isotopes have nearly constant 
E(2+) (~ 1200 keV) values for A=102-130. This is at- 
tributed to the near degeneracy of the neutron lgj/ 2 and 
2d 5 / 2 single particle orbitals which enhances pairing cor- 
relations. Similarly, there is a general belief that for Sn 
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isotopes beyond 132 Sn, there will be enhancement of pair- 
ing resulting in a constant 0^ — 2^ spacing, as revealed 
in calculations with realistic interactions [HI [13] ■ How- 
ever, suddenly at N=84, for 134 Sn, the E(2f) decreases 
to 726 keV. This decrease of £(2+) (Fig. g]) is clearly 
correlated to the decrease in pairing manifested through 
the reduction in A„ (~ 1.2 MeV at N=81 to 0.6 MeV at 
N=83). This observation is strengthened by the most re- 
cent precise mass measurements [18|]. On the other hand, 
for Pb isotopes above N=126, the A„ values are very 
similar to those below N=126 (~ 0.68,0.32,0.62 MeV at 
N=123,125, 127, respectively) and so a near constancy of 
the £(2+)s below and above N=126 closure (0.899, 0.803, 
0.800 MeV for N=122, 124, 128, respectively) has been 
observed. Pairing for different isotopes of Te, Xe and 
Ba nuclei has been estimated. At N=82, E(2f) (1.279, 
1.313 and 1.435 MeV, respectively) for Te, Xe and Ba 
is ~ 2A P (A p ~ 0.66, 0.82 and 0.87 MeV, respectively). 
Compared to that in N=78 (A„ ~ 1.1 MeV for Te, Xe 
and Ba), beyond N=82, reduction in neutron pairing (A„ 
~ 0.72, 0.83, 0.84 MeV at N=84 for Te, Xe and Ba, re- 
spectively) is observed. It gives rise to depressed E(2^) 
(0.606, 0.589 and 0.602 MeV, respectively) and relatively 
hindered B(E2)s depending on the degree of reduction. 
Weaker A p at N=82 favour easy generation of 2^~ state 
without any need for neutron shell erosion. 

Analysing the empirical observations of the features of 
neutron-rich nuclei one can frame the following rules. 

• For proton - magic nuclei, reduction in neutron 
pairing will give rise to 

— anomalous decrease in E(2^~) with slower in- 
crease in deformation for isotopes with neu- 
trons in unfilled subshells. For normal pairing 
observed near stability, similar value of low 
E(2+)'s will indicate higher value of B(E2) or 
larger deformation. 

— appearance of new shell closures manifested 
through enhancement of E(2+) for isotopes 
with filled up sub-shells. The energy gap be- 
tween the sub - shells must be substantially 
greater than the neutron pairing energy. 

• For nuclei with a few valence particles near shell 
closure having unfilled proton sub-shell, reduced 
neutron pairing will lead to 

— anomalous decrease in E(2^~) with slower in- 
crease in deformation for isotopes with neu- 
trons in unfilled sub-shells, 

— onset of deformation for N mQg i C nucleus, re- 
sulting in erosion of shell gap if proton pair- 
ing shows strong enhancement. This will be 
usually observed for lighter mass nuclei. 

To have a quantitative estimation of the effect of pair- 
ing on depression in the E(2~£) values and on the onset 
of deformation, shell model calculation has been done for 
oxygen isotopes. It has been shown that with reduction 
of + diagonal two body matrix elements of Bonn A real- 
istic interaction [l9j], E(2f) energies of 18 ~ 22 decreases 
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FIG. 5: (a) The equi- E(2+) energy contours are shown as 
a function of deformation (/?) and neutron pairing gap (An), 
(b) The variations in E(2+ ) energies for 134-140 g n with re- 
duction in neutron pairing gap. 



without any increase in the B(E2) values with improve- 
ment in prediction of a shell gap at N=14, usually not 
predicted by Bonn-A. The same is also true for even-even 
Sn isotopes beyond 132 Sn as also discussed in Ref. (l6| . 

From cranked Hartree-Fock-Bogoliubov (CHFB) cal- 
culations with a pairing plus quadrupole Hamiltonian, 
equi-energy contours for E(2f) energies in 124 JTe, well 
studied within this formalism [2(| have been plotted 
in Fig. [5^i, for variations of neutron pairing gap and 
quadrupole deformation. For all these contours the pro- 
ton pairing gap has been kept constant at 1.5 MeV, which 
is obtained from the experimental odd-even mass differ- 
ence. These contours are useful to understand the effect 
of variation of pairing on E(2f ) for fixed value of defor- 
mation and vice- versa. For example, the measured F,(2^) 
value (~ 350 keV) for 124 Xe can correspond to two ex- 
treme combinations of neutron pairing gap and deforma- 
tion. They are (A„ = 1.4 and ,3=0.30) and (A„=0.5 and 
,3=0.11) (Fig. [Sk). It implies that same depressed value 
of energy can result either due to increase in deformation 
with usual pairing or relatively smaller deformation with 
reduced of pairing. 

Self-consistent solutions have been obtained to inves- 
tigate how deformation changes with decreasing pair- 
ing and vice - versa for 12i Xe. The results are shown 
in Table I. It is found that even without any change 
of quadrupole interaction strength, reduction in neu- 
tron pairing strength favours a slow onset of deforma- 
tion (Part (a) of Table I). It is also evident that same 
depressed E(2^~) energy (e.g., ~ 120-150 keV) can re- 
sult either from a reduced pairing (A„ ~ 0) and rela- 
tively smaller deformation (/?2 — 0.26) (Table la), or from 
enhanced deformation (/?2 — 0.37) for stronger pairing 
(A„ ~ 0.6 MeV) (Table lb). 

To test the observation for proton-magic nuclei, we 
have done another set of calculations for Sn isotopes 
above s Sn. Single particle energies have been obtained 
from shell model Hamiltonian in Refs. [8j, [lg, l2l| . The 
single particle energies for the wlgg/2 and vlhn/2 orbitals 
have been chosen with sufficient care so that Sn isotopes 
above A=132 do not show any protons or neutron exci- 
tations across the Z=50 and N=82 shell closures, respec- 
tively. The remaining single particle energies are calcu- 
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TABLE I: Effects of (a) decrease in neutron pairing strength 
(gn) and (b) increase in quadrupole interaction strength, X2, 
on ground state deformation, pairing gaps and E(2f) values 
of 124 Xe. In part (a), the X2 value is kept constant at 72.0, in 
(b), the value of g n is 20.0. The proton pairing strength (g p ) 
is kept constant at 26.0 in all these calculations. The energies 
are in MeV. 



g n 


fa 


A p 


A„ 


E(2+) 


(a) 25 


0.087 


1.368 


1.875 


0.773 


20 


0.182 


1.222 


1.047 


0.304 


15 


0.230 


1.171 


0.433 


0.172 


10 


0.261 


1.168 


0.00003 


0.125 


5 


0.261 


1.168 


0.00001 


0.114 


X2 


fa 


A p 


An 


E(2+) 


(b) 72 


0.182 


1.222 


1.047 


0.304 


77 


0.256 


1.167 


0.834 


0.233 


82 


0.366 


1.116 


0.605 


0.151 



lated using Nilsson prescription |22j | with proper normal- 
isation. 

From the Fig. [5}d, it is observed that as pairing energy 
for neutron decreases, the shell closure at N=90 |21j man- 
ifested through sudden jump in E(2^~) energy of 140 Sn 
compared to that in 138 Sn becomes more pronounced. 
For stronger pairing the E(2j*~) energies of semi-magic 
nuclei show much slower variation with increasing neu- 
tron numbers, as also observed for isotopes of Sn below 



132 Sn. 

From the experimental data on odd-even staggering of 
masses and SM and CHFB calculations, we have shown 
that pairing plays an important role in the shell evo- 
lution of neutron-rich nuclei. We have identified a few 
distinctive features of neutron-rich nuclei based on our 
analysis of the empirical data. For proton-magic nuclei, 
due to reduction of neutron pairing away from stability, 
each sub-shell closure will be manifested as a new shell 
closure, if energy gap between sub-shells is greater than 
A„. For other neutron- rich isotopes for which proton 
numbers are not magic, the E(2^~) energies at each ex- 
isting neutron shell closure will show dramatic decrease 
depending on the extent of reduction in neutron pairing. 
This reduction will also be manifested through slow onset 
of deformation at shell closures. But these values of the 
deformation will not correspond to the same variation as 
seen on the stability. The increase in deformation with 
decrease in energy of the 2^ state will be much slower 
for reduced pairing. For lighter mass nuclei, increase in 
proton pairing has been found to play a significant role 
in the observed shell erosion. This work therefore conclu- 
sively shows the role of variation in pairing as a function 
N in nuclear shell evolution of neutron - rich nuclei. This 
method of analysis is also applicable to exotic proton-rich 
nuclei. The microscopic origin of variation in pairing with 
N constitutes a separate important issue. Our work also 
warrants more precise measurement of masses away from 
stability and analysis of mass data to search for newer 
isotopes for which such shell evolution can be expected. 
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